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I review the main characteristics of structure formation in the quintessen¬ 
tial Universe. Assuming equation of state w = p/g =const I provide a brief 
description of the background cosmology and discuss the linear growth of 
density perturbations, the strongly nonlinear evolution, the power spectra 
and rms fluctuations as well as mass functions focusing on the three values 
w = — 1, — 2/3 and —1/3. Finally I describe the presently available and 
future constraints on w. 

PACS numbers: 95.35.+d, 98.62.-g, 98.62.Ai, 98.62.Ck, 98.80.-k, 98.80.Es 


1. Introduction 

Our knowledge of background cosmology has recently improved dra¬ 
matically due to new supernovae and cosmic microwave background data. 
Current observations favor a flat Universe with matter density flo = 0.3 [1] 
and the remaining contribution in the form of cosmological constant [2, 3] 
or some other form of dark energy. The models with cosmological constant 
are known, however, to suffer from two major problems. One is related to 
the origin of the constant - it cannot be explained in terms of the vacuum 
energy since its energy is orders of magnitude smaller. The other is the 
lack of explanation why the present densities in matter and cosmological 
constant are comparable. 

A new class of models that solve these problems and also satisfy present 
observational constraints has been proposed a few years ago [4]. In these 
models the cosmological constant is replaced with a new energy component, 
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called quintessence, characterized by the equation of state p/g = w / — 1 . 
The component can cluster on largest scales and therefore affect the mass 
power spectrum [5] and microwave background anisotropies [6, 7]. 

The investigations of the physical basis for the existence of such com¬ 
ponent are now more than a decade old [8]. One of the promising models 
is based on so-called “tracker fields” that display an attractor-like behavior 
causing the energy density of quintessence to follow the radiation density in 
the radiation dominated era but dominate over matter density after matter- 
radiation equality [9, 10]. It is still debated, however, how w should depend 
on time, and whether its redshift dependence can be reliably determined 
observationally [11, 12, 13]. 

A considerable effort has gone into attempts to put constraints on models 
with quintessence and presently the values of — 1 < w < —0.6 seem most 
feasible observationally [14, 15]. Here I review the main characteristics 
of structure formation in the quintessential Universe which may provide 
constraints on the equation of state. In the last Section I discuss the current 
status of observational limits on w and future perspectives. 

2. Background cosmology 

Quintessence obeys the following equation of state relating its density 
qq and pressure pq 


PQ = wqq, where — 1 < w < 0. 


( 1 ) 


The case of w = — 1 corresponds to the usually defined cosmological con¬ 
stant. 

The evolution of the scale factor a = R/Rq = 1/(1 + z) (normalized to 
unity at present, z is the redshift) in the quintessential Universe is governed 
by the Friedmann equation 


da Hq 
d t u(a ) 


( 2 ) 


where 


u(a) 


i + n 0 



+ qo 


l 


,l+3ui 



( 3 ) 


and Hq is the present value of the Hubble parameter. The quantities with 
subscript 0 here and below denote the present values. The parameter 11 is 
the standard measure of the amount of matter in units of critical density 
and q measures the density of quintessence in the same units 



£?crit 


(4) 
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Fig. 1. Left panel: Time evolution of the scale factor in different models. Right 
panel: The present age of Universe in units of H^ 1 in flat models as a function of 
Uo for different w. 


The Einstein equation for acceleration d 2 a/dt 2 = —4irGa(p + g/ 3) shows 
that w < —1/3 is needed for the accelerated expansion to occur. The left 
panel of Fig. 1 shows the evolution of scale factor in different models. The 
right panel presents the dependence on w of the present age of Universe 


to 


1 A 1 

— / u{a)da. 
Flo JO 


( 5 ) 


Solving the equation for the conservation of energy d(gQa?) / da = —3 pQa 2 
with condition (1) we get the evolution of the density of quintessence which 
for w = const, the case considered in this paper, reduces to 


GQ 


GQ, 0 a 


3(l+ui) 


( 6 ) 


The evolution of and q with scale factor is given by 


= SW(a) 

a 


q(a) 


qoU 2 {a ) 
a l+3w 


( 7 ) 


while the Hubble parameter itself evolves so that H(a ) = Hq/[ a u(a)]. 


3. Linear growth of perturbations 

The linear evolution of the matter density contrast 5 = Sg/ g is governed 
by equation [16] 

5 + 2—5 — 47t GgS = 0 
a 


(8) 
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Fig. 2. The linear growth rate of density fluctuations for fio = 0.3, qo = 0.7 in 
three cases of w = —1, —2/3 and —1/3. 


where dots represent derivatives with respect to time. For flat models and 
arbitrary w an analytical expression for D(a), the growing mode of the 
time-dependent part of 5, was found [17]. With our notation and the nor¬ 
malization of D(a) = a for Q = 1 and q = 0 it becomes 


D(a) = a 2 Fi 


'_L w ~ 1 i _ A j-3^ 1 ~ 

3 w ’ 2iv ’ 6w ’ flo 


( 9 ) 


where 2-^1 is a hypergeometric function. The solutions (9) for different w 
and cosmological parameters Qo = 0.3 and qo = 0.7 are plotted in Fig. 2. 
The peculiar velocity field in linear perturbation theory is obtained from 


2 / 

3 QH g 


( 10 ) 


where g = —'V(j)/a is the peculiar gravitational acceleration and / is the 
dimensionless velocity factor 


a D 
a D 


( 11 ) 


For flat models this formula can be evaluated analytically using Eq. (9). The 
dependence of / on Qq at present (z = 0) for flat models with different w is 
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Fig. 3. Left panel: The velocity factor / at present as a function of fio for flat 
models with different w. Right panel: The redshift dependence of / for Qq = 0.3, 
<7o = 0.7 and different w. 

shown in the left panel of Fig. 3. We see immediately that the dependence on 
w is very weak. However, as shown in the right panel of Fig. 3, when going 
to higher redshifts we find that the velocity factor is much more sensitive to 
w which gives some hope for applying it to determine w from local peculiar 
velocity held. 

4. Strongly nonlinear evolution 

The simplest model of formation of bound objects (called the spherical or 
top hat model) [18, 19, 20, 21] describes the nonlinear evolution of a uniform 
spherical density perturbation in an otherwise homogeneous Universe. If the 
overdensity is big enough the perturbed region will slow down the expansion 
and eventually turn around and collapse. The evolution of the size of such 
a perturbation can be described [22] by the Einstein equation 

r 1 \ 11 

- = —4t tG [w + -j qq + -Q ■ (12) 

Together with the Friedmann equation (2) it can be solved to obtain the 
ratio of the cluster density to the background density at the time of turn¬ 
around 

C = = (^) (13) 

Qb V 4 J 

where the last formula is the approximation given in [22]. 
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fig 2-con 


Fig. 4. Left panel: dependence of A c on flo and w for flat models and collapse time 
Zcoii = 0. Right panel: A c as a function of z co n for flo = 0-3, qo = 0.7 and different 

w. 


One of the most useful quantities one can derive from the spherical top 
hat model is the density of the virialized halo in units of the critical density 
at the time of collapse corresponding to redshift z co \\ 

Ac = (Zcoll) = C(*a)fi(3coll) f— V • (14) 

£>crit \r co \\J V 1 + ^coll / 

The redshifts z co ii and zt a are related by the assumption that the collapse 
time tcolb the time corresponding to r —> 0, is twice the turn-around time 

1 rata 

ha = tt / u(a)da. (15) 

ti o Jo 

In reality the object does not decrease its size to zero but virializes to reach 
the effective final radius r co u. The ratio by which the halo collapses is 
estimated from the virial theorem which gives [22] 

—— = 1 ~ ^ /2 , , (16) 

rta 2 + Tjt- 3 ? 7„/2 

where 


2 q(zta) 

C hl(z u ) 

2 g(^coll) / 1 + ZcoiA 3 

C 12(-2-coll) \ 1 T -2-ta ) 


(17) 

(18) 


Figure 4 shows how A c depends on , w and the redshift of collapse. 
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5. Power spectrum of density fluctuations 

Power spectrum P(k, a) is defined as the Fourier transform of the cor¬ 
relation function of density fluctuations 

P(k,a ) = J £{r,a) e _lkr d'V. (19) 

The spectra for Universe dominated by cold dark matter (CDM) have been 
widely discussed in the literature, e.g. [23]. For the present time (a = 1) 
the power spectrum is usually written in the form 

P(k) = Ak n T\k ) (20) 

where n measures the slope of the primordial power spectrum (we will as¬ 
sume n = 1), T is the transfer function and A is a normalization constant. 
In the presence of cosmological constant (ACDM) the transfer function T\ 
can be approximated by [23] 

T l( p ) = hl 34p)2 4 ^ [1 + 3.89p+ (16.Ip) 2 + (5.46p) 3 + (6.71p) 4 ]~ 1/2 , (21) 

where p = /c/(r/iMpc _1 ), T = flofaexp[—Db(l + V2h/klo)\ and h = #o/[100 
km/(s Mpc)] = 0.7. We assume the barion contribution to the density 
parameter fib = 0.04. 

For flat models with quintessence the modification of the transfer func¬ 
tion has been proposed by Ma et al. [5]. The transfer function in equation 
(20) is then Tq = Tq\T\, where Tq\ = Tq/T\ is approximated by fits 
given in [5]. The present linear power spectra obtained for w = — 1,— 2/3 
and —1/3 assuming COBE normalization are shown in the left panel of 
Fig. 5. 

The rms density fluctuation, <r, at comoving smoothing scale R is given 
by 

^ = ( 2 bfj d3kP ( k ) W T H (kR) ( 22 ) 

where the smoothing is performed with the top hat filter 

^™ (M) = w(^~ C 0 SWi )' (23) 

The dependence of a on smoothing scale for flat models with different w 
is shown in the right panel of Fig. 5. A particularly useful quantity, con¬ 
strained by cluster abundance is the rms fluctuation at the scale of 8/i -1 
Mpc. Due to COBE normalization its values turn out to depend strongly 
on w and we get as = 0.96,0.80 and 0.46 for w = —1,—2/3 and —1/3 
respectively. 
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Fig. 5. The linear power spectra (left panel) and rms fluctuation as a function of 
smoothing scale (right panel) for flat models with different w. 


6. Mass functions 

One of the most important measures of structure formation is provided 
by the mass function of bound objects. Using the analytical prescription 
of Press and Schechter [24], we can estimate the cumulative mass function 
(the comoving number density of objects of mass grater than M ) 

pO O 

N(> M ) = / n(M)dM (24) 

Jm 


where n(M ) is the number density of objects with mass between M and 


M + d M 


n(M ) 


2\ 1/2 gb S c da /_ 

7 t J M a 2 d M ^ y 2a 2 J 


(25) 


In the expression above, g b is the background density, a is the rms density 
fluctuation at comoving smoothing scale R described in Section 5. The mass 
is related to the smoothing scale by M = 4Trgy ) R 3 /3. The parameter 5 C is 
the halo density at collapse as predicted by linear theory. We adopt here its 
standard value 5 C = 1.69 (exact for the Einstein-de Sitter Universe) since it 
depends very weakly on Uo and w [22], 

Figure 6 shows the cumulative mass functions calculated from equation 
(25) with flo = 0.3 and qo = 0.7 for three models with w = —1, —2/3 and 
— 1/3. The strong dependence on w is due to the COBE normalization of 
the power spectra. The Figure also shows data for rich clusters of galaxies 
from [25]. For our choice of cosmological parameters the value of w ~ — 1 is 
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Fig. 6. The Press-Schechter cumulative mass functions for different w assuming 
O 0 = 0.3 and q 0 = 0.7. 


preferred. However, the Press-Schechter formulae are known to underesti¬ 
mate the mass functions on the scale of clusters of galaxies when compared 
to N-body simulations. The mass functions are also very sensitive to the 
value of Ho- Therefore when more exact predictions are used and Ho is 
measured to be close to 0.3 with a small error bar, we can expect higher 
values of w to be preferred by the cluster data. 

7. Constraints on w 

The primary constraint on the value of w comes from the observations 
of accelerating expansion, which can only be obtained in models with w < 
— 1/3. Age of Universe is quite sensitive to w and increases for lower w, 
however the accuracy of our knowledge of both Hq and to is not good enough 
to put strong constraints on w. Current estimates are consistent with w < 
-0.5. 

One of the strongest arguments for the existence of dark energy comes 
from the studies of cosmic microwave background (CMB). Although its 
power spectrum is weakly sensitive to w (e.g. the height of the first acoustic 
peak is somewhat increased and its position is shifted to higher multipoles 
for lower w), it has already provided some limits. Combining the data from 
COBE and recent balloon experiments Balbi et al. [6] find —1 < w < —0.6 
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while Baccigalupi et al. [7] estimate the best-fitting value of w to be w = 
—0.8. The ongoing and future satellite experiments are expected to put 
even stronger limits on w. 

Among the promising probes of dark energy are also supernovae la. The 
comoving distance they serve to measure is sensitive only to the interesting 
cosmological parameters and the errors related to supernova evolution or 
extinction are estimated to be small. The existing data restrict w only 
weakly [26], but future experiments like The Supernova Acceleration Probe 
are expected to measure the value of w with a few percent accuracy [27]. 

Structure formation also offers methods to constrain the cosmic equation 
of state. The suppression of linear growth of density fluctuations for higher 
w alone shows that only for w < —1/2 the structure observed today could 
have evolved from small initial perturbations deduced from CMB observa¬ 
tions. The same range of acceptable w values follows from the behaviour of 
<7g which is a strongly decreasing function of w. 

The most promising tests are based on the number counts of galaxy 
clusters. It turns out [22] that the slope of comoving abundance as a function 
of redshift depends sensitively on w and therefore can be used to break 
degeneracies between w and other cosmological parameters that appear e.g. 
in the analysis of CMB. Such measurements are expected to be performed 
using the proposed new X-ray and Sunyaev-Zeldovich effect surveys [28] and 
the ongoing DEEP Redshift Survey [29]. The constraints from structure 
formation appear to be complementary to those from supernovae and CMB 
measurements. 
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